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Abstract 
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has been used to prove that: i) the quantum bounce that replaces the big bang singularity is generic; 
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times had to be semiclassical before the bounce. Here we consider a family of exact solutions 
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considerations. These turn out to be enough to select, from all the generalized coherent states, 
those that behave semiclassical at late times. We study then the properties of such states near 
the bounce where the most 'quantum behavior' is expected. As it turns out, the states remain 
sharply peaked and semiclassical at the bounce and the dynamics is very well approximated by 
the 'effective theory' throughout the time evolution. We compare the semiclassicality properties 
of squeezed states to those of the Gaussian semiclassical states and conclude that the Gaussians 
are better behaved. In particular, the asymmetry in the relative fluctuations before and after the 
bounce are negligible, thus ruling out claims of so called 'cosmic forgetfulness'. 
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I. INTRODUCTION 



Loop quantum cosmology has provided in past years a useful framework to ask questions 
about the quantum nature of the early universe [U [2]. Closely related to loop quantum 
gravity [3], the formalism has shown that, for isotropic models, the big bang singularity is 
replaced by a quantum bounce |1H6]. These results have been obtained for closed and open 
FRW cosmologies [7[ E] with and without a cosmological constant [10] for a massless scalar 
field, and recently extended to the massive case |11]. Recently, some of these results have 
been extended to anisotropic [T2HT5] and some inhomogeneous models [16]. In all these 
cases, quantum gravitational effects have been shown to exert a repulsive force and halt 
the collapse and launch an expanding superinflationary phase before the quantum gravity 
effects die out and the universe then follows the standard general relativity (GR) dynamics. 

An exactly solvable model for the k=0 FRW with a massless scalar field has allowed to 
prove analytically the generic nature of the quantum bounce [T7]. In particular, it was shown 
that the bounce occurs for all physical states, and the energy density possesses a supremum 
for physical states [H] . Regarding the fluctuations, it has been shown that they are bounded 
across the bounce implying that a semiclassical state at late times after the bounce had to 
come from a semiclassical state before the bounce [18]. Recently, Kaminski and Pawlowski 
have imposed stronger bounds on the asymmetry in relative dispersions before and after the 
bounce [T9] , thus very strongly ruling out some claims of loss of semiclassicality in [20l |2T] . 
These robustness result have support on the fact that, within a large possibility of 'loop 
quantizations', physical considerations select a unique consistent quantization [22] (the one 
introduced in [B]). 

An interesting feature of this system is that one can approximate very well the dynam- 
ics followed by semiclassical states by an effective description [6j. That is, the classical 
limit of the dynamics is not that of GR + scalar field but rather is defined by an effective 
Hamiltonian. There have been several approaches to reach this result. For instance, in [23] 
the author considered kinematical coherent states and used them to compute the effective 
Hamiltonian as the expectation value of the quantum Hamiltonian constraint of LQC. More 
recently, the authors of [24j showed that, in the path integral description of the model, the 
paths that contribute the most are those that satisfy the effective equations and not those 
that follow the classical GR equations. 

The purpose of this paper is to construct exact physical semiclassical states for the solv- 
able model [17j. The strategy will be to define a class of coherent initial states and, using the 
analytical control we have on the space of solutions to the quantum constraint equations, 
produce exact physical coherent states, much in the spirit of [25]. Then, we impose well 
motivated conditions on the behavior of the fluctuations of physical operators to select from 
all possible coherent states, those that exhibit a semiclassical behavior. It turns out that 
these conditions are sufficient to fix enough parameters of the states and to have states that 
are 'peaked' on a given (physical) phase space point. Having found semiclassical states, we 
explore the behavior of these states near the bounce, where one expects the state to be most 
quantum. In particular we consider the relative fluctuations of volume and find that, for all 
states and for all times, the relative fluctuations are bounded by their asymptotic value. In 
fact, the minimum value turns out to be precisely at the bounce. The expectation value for 
volume and density on semiclassical states are very well approximated by the values given 
by the effective theory. This leads us to conclude that the effective description is indeed a 
very good approximation to semiclassical physical states. A natural question is whether the 



2 



class of states considered is the most general one. For that purpose, we generalize the class 
of Gaussian states to the so called squeezed states. The idea is to explore the parameter 
space near the Gaussian states -that are well behaved semiclassical states- and see whether 
these squeezed states exhibit better semiclassical properties. 

Of particular interest in this regard is the asymmetry in volume fluctuations. Claims of a 
loss of coherence, or 'cosmic forgetfulness' arise due to results in an approximate framework 
^20j (motivated by, but not within LQC) that suggest that the fluctuations across the bounce 
might change significantly and spoil the semiclassical nature of the state [201 [2lj . Even when 
strong bounds against this possibility within LQC have already appeared [19], our exact 
formalism is particularly well suited for addressing this question. By imposing bounds on the 
total asymptotic dispersion in volume (at late or early times with respect to the bounce), we 
are able to bound how far we can be from Gaussian states and still have a semiclassical state. 
The range on parameter space allowed turns out to depend on the physical requirements for 
the state, in such a way that, for 'large volume universes', it is severely constrained. This is 
turn sets very strong bounds on the allowed asymmetry in volume fluctuations, consistent 
with [191 ^^^y in particular, disproving some of the statements of [201 EI] that claim the loss 
of coherence across the bounce for semiclassical states. 

The structure of the paper is as follows. In Sec. |TT] we recall the solvable model in loop 
quantum cosmology together with the basic operators to be studied. In Sec. Ill we introduce 
the Gaussian states and compute expectation values and fluctuations for the physical oper- 
ators. Semiclassicality conditions are explored in Sec. IV, together with the behavior of the 
states near the bounce. Squeezed states are introduced and compared to the Gaussian states 
in Sec. |V} We end with a discussion in Sec VI There are four Appendices. In Appendix [A] 
we study the symplectic reduction of the 'effective theory' to isolate the physical degrees of 
freedom. In Appendix |B] we recall some useful integrals and formulas. In Appendix [C] we 
present an analysis of the errors introduced by considering Gaussian states, showing that 
they are severely suppressed and justify our analysis in the main body. In Appendix |D] we 
present in a table a comparison of different quantities for two families of Gaussian states. 
Some of the results reported here were summarized in the short manuscript [26J, so it is 
natural that there is some overlap of part of the material. 



II. SOLVABLE LOOP QUANTUM COSMOLOGY (SLQC) 

In this section we shall recall the model we are considering, namely, the solvable k=0 FRW 
model coupled to a massless scalar field. Our presentation will be self-contained and slightly 
different from that of [T7]. In terms of the phase space variables used in loop quantum 
gravity (LQG) [3], namely a connection and a densitized triad E"", the homogeneous 
gravitational sector can be expressed as, 

where (°Co'*,°e") are a set of orthonormal co-triads and triads compatible with the fiducial 
(flat for k = 0) metric °qab and Vo is the volume of the fiducial cell, introduced to define the 
symplectic structure, with respect to °qab- The phase space is characterized by conjugate 
variables {c,p) satisfying {c,p} = (87rG7)/3. where 7 ~ 0.2375 is the Barbero-Immirzi 
parameter used in LQC. The triad p is related to the physical volume V of the fiducial cell, 
as \p\ = = Vo^^a'^ where a is the scale factor and, on the space of classical solutions. 
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c = 70. It is convenient to introduce new variables |17] . 

/3:=ec/p^/2 and = £p^/V(27r4i7) (2.2) 

where h{P, u} = 2^. Here e = ±1 is the orientation of the triad with respect to that of "w*. 
The classical constraint, with the choice N = V, then becomes 

where the matter content consists of a massless scalar field with canonical momenta 
Ptp. Thus, the kinematical phase space can be though of as a four dimensional space with 
coordinates {/3,V, (j),p^). The strategy for quantization is to define an operator associated 
to C and to look for states that are annihilated by the constraint operator. 
In the (/3, 0) representation, the operators basic are represented as 

/3 ■ 0) = /3 Xi(3, ^) ; V- Xi(3, 4>) = -^Wpi^ 0) ; (2.4) 
P^-xiPA) = -ih^xiPA)- (2.5) 

However, the operator (3 is not well defined since /3 is periodic with period 27r/A, so it has to 

be replaced by a well defined operator (3 := sin(A/3)/A (see [T7j for details). The quantum 
constraint of sLQC becomes then [T7] . 

■x{M)=a^{'^l\-x{M) (2.6) 



-^^"'^^ \ X dp 

with (5 G (0,7r/A), a := \/12t:G. In the loop quantum cosmology literature, the value of A 
is chosen such that A^ = A corresponds to the minimum eigenvalue of the area operator in 
loop quantum gravity (corresponding to an edge of 'spin 1/2'). With this choice the free 
parameter A becomes A^ = 4-\/37r7£pi. By introducing 

X = a-Mn|tan(A/3/2)| (2.7) 

the quantum constraint can be rewritten in a Klein-Gordon form 

dlx{xA) = dlx{xA) ■ (2.8) 

As usual for this system, in its classical evolution the scalar field is a monotonic function 
and can play the role of internal clock. In the quantum theory, one can also think of the 
evolution of the state with respect to 0. A general solution x(x,0) to Eq. (2.8) can be 
decomposed in the left and right moving components: 

X = X+(0 + 2;) + X-(0 - x) := x+i^+) + X-i^^) . (2.9) 



The physical states that we shall consider are positive frequency solutions of (2.8). Since 
there are no fermions in the model, the orientations of the triad are indistinguishable and 



Here we have adopted the notation for the variable /3 introduced in [55] and used thereafter. 
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x{x, 0) satisfy the symmetry requirement x{~^: 0) = ~x{^j 0)- Thus, we can write 0) = 
- F(x_))/\/2, where F is an arbitrary 'positive frequency solution'^. The physical 
inner product on solutions is given as [T7] 



phy 



[Xiix,(f))d^X2{x,(f)) - ((9<^xi(a;,0))X2(a;,0)]dx (2.10) 

/oo 

[9,.Fi(x+)F2(x+) - 9,.Fi(x.)F2(x_)] dx . (2.11) 
-oo 



We can now compute the expectation values and fluctuations of fundamental operator such 
as V\^^,x and p^, where x is related to the operator (3. For any state on the physical Hilbert 
space the expectation value of the volume operator at 'time 0' is given by 



{y)<^ ■= ix, y\<pX)phj = 27r74i(x, l'>lx)phy 
where |z>| is the absolute value operator obtained from 

u = cosh.(ax)idx ■ 

a 



(2.12) 



(2.13) 



Using the inner product (2.11) the expectation value of |z/| is given by 

/oo 
[9,F(x+)(z>F(x+)) - 9,F(x_)(-z>F(x_))] dx 
-oo 

o\ r<x> 

[(9^F(x+) cosh{ax)dxF{x+) + d^F^x^) cosh.{ax)drcF (x^)] dx 



2A 

a 
4A 

a 



dF 



dx 



cosh(a(x — (p)) dx . 



(2.14) 



Is easy to check that |z>| = |i>|"''. From these expressions one can find the expectation value 
of V'^ and the dispersion of the operator, 



{AV)l = Yo + Y+ e^"^ + r_ e-2"<^ , 
with V±, Wq, W±, Yq and Y± being real and positive, given by 



(2.15) 
(2.16) 
(2.17) 



Wo 



47r74iA 



a 



dF 



dx 



dx 



3G 

„2M \2 



d^F dF dF d^F 



dx 



3G ' 



dx^ dx dx dx^ 

d^F dF dF d^F 



dx^ dx dx dx^ 



dx 



^ To be precise, F{x) is a positive momentum function, i.e. with a Fourier transform that has support on 
the positive axis. With such a choice, the solution to the constraint equation become of positive frequency. 
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■VI 



for normalized states. Here Yq = Wq — 2V+V- and Y± = 
for Wq an d W± correct those found in |18J. 

From (2.15^, it follows that the expectation value of the volume V 



Note that the expressions 



is large at both 
2{V+V_ 



^1/2, 



\X\ 



very early and late times and has a non-zero global minimum V^in 

The bounce ocurrs at time 0^ = (2a)^^ \n(y^/V+) [17]. Around = 0^, the expectation 
value of the volume {V)(f, is symmetric. Similarly, {v)(j> is symmetric across the value 
4>^^ = {4a)~^\n{W-/W^) for the scalar field. A trivial observation is that if 0^ = 0^^, 
difference in the asymptotic values of the relative fluctuation 



the 



D :-- 



lim 



(AV) 




w. w. 



VI 



(2.18) 



vanishes. It should be noted that this quantity quantifies the change in semiclassicality 
across the bounce as pointed out in [T8] . 

The other observable that one might want to consider as fundamental, is an observable 
naturally conjugate to p^. Given that we are considering physical states, this does not 
correspond to the operator as happens in the kinematical setting. Instead, as illustrated 
in the Appendix |X| the quantity that is conjugate to p^, in the reduced physical model, 
corresponds to the quantity x(/3). Thus, it is natural to consider the family of operators x^, 
corresponding to the observable 'x at time 0'. The expectation value of x is defined as. 



(2.19) 



(Xl,^^X2)phy = 2[(l*IXl,X2)phy + (Xl, |^^|X2)phy] 

where |x|F(x+) = xF(x+) and |x|F(x_) = — xF(x_). We have included the two terms in 
the definition (2.19) because |x| is not symmetric. We can now compute the expectation 



values and fluctuations of x for any state of the physical Hilbert space 





= X, 


-0, 




= X2 


- 20Xi + 02 


(Ax) 


= X2 





with Xi and X2 real (state dependent) constants, given by 

/■oo 

Xi = -2i I F{x)xd^F{x)dx - 

F{x)x d-^F{x)dx-2i I |F(x)pxdx', 



Xo 



-2i 



-00 
00 



|F(x)rdx 



(2.20) 
(2.21) 
(2.22) 



(2.23) 
(2.24) 



for normalized states. Is important to note that (x)<^ oc 0, (x^),^ oc 0^ but the fluctuation 
Ax is independent of 0! This property indeed gives support to the proposal that x be the 
conjugate of p^. Finally, we compute the expectation values for p^, a Dirac observable as 



-ihdrh 



{P4>) 
(Pi) 



2h 



dF 



dx 



-2h^i 



dx 

dF d^F 
dx dx^ 



dx', 



(2.25) 
(2.26) 
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for normalized states. In the following sections, it will be convenient to consider a general 
physical state written as, 



F(A:) e-''=('^+") dA; - / e-''^^^"") dfc 



(2.27) 



where the Fourier transf orm F{k) contains all the information of the state. Positive fre- 
quency solutions to Eq. (2.8) means that F{k) has support on positive fc's only. Further- 
more, we can consider initial states, for an 'initial time' = 0, in order to compute some 
of the quantities that determine the expectation value of relevant operators. It then suffices 
to specify F{k) and from there to compute all the quantities V±,Wo,W±, Xi, X2, (Pfj,) and 



III. GAUSSIAN INITIAL STATES 



In this section we shall consider physical states with initial states given generalized Gaus- 
sian states, defined by functions F{k) of the form: 



Fik) 



J^n ^~ (k- ko)'^ / a'^ ^ipk 

0, 



for k > 0, 
for < 0, 



(3.1) 



with cr > 0, fco > 0, n = 0, 1, 2, .... That is, we are choosing Gaussian states centered around 
the point ko, with 'dispersion' given by a. 

The physical states one has to consider are vanishing in the negative k axis, to have 
positive frequency solutions. In order to gain analytical control over all the quantities at 
hand, we shall in what follows consider instead the Gaussian states for all values of k. Thus, 
the quantities we shall compute are an approximation to the real quantities but, as we shall 
argue in detail, this approximation is justified when the values of k^ and a are such that 
the wave function has a negligible contribution from the negative axis. A detailed analysis 
of the errors introduced by this simplification can be found in the Appendix [Cj Therefore, 
from now on, in order to have states that approximate positive frequency solutions, we shall 
impose the condition ko ^ a. Further consistency conditions on a, motivated by physical 
considerations will be derived below. 

The norm of the state in k space is 

/•oo _ roo 

(X,x)phy = -2i / F{k)ikF{k)dk = 2 / k\F{k)\^<lk. 



If we take the initial states in the fc-space as in (3.1), the norm is, 

(X,x)phy = 2/ ke-e-'^'-'^^'I'^'dk. 



which can be approximated by (x, x)phy = 2 kk'^'^e '^^^ ^"^^/'^^dA; (see Appendix 
justification of this approximation). 

With the change of variables u = ^{k — ko) this last integral takes the form 

2n+l 



c 



for a 



Ixllphy 2 



cr 


2n+2 


POO 


V2ko 

u H 

a 


[V2\ 




' —00 



e"" du . 



(3.2) 
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For concreteness, we shall consider the state corresponding to n = 0. In this case 

e-"'dM = v^A;oa. (3.3) 

The normalized n=0 states can be used to compute explicitly the expectations values 
and fluctuations of several physically interesting operators. The corresponding results for 
n=l are summarized in the Appendix [D| 



Ixllphy 



a 


2 /•oo 


[V2\ 


J — oo 



u + 



a 



A. Expectation Values for Basic Observables 



Let us explore the basic observables for the Gaussian states in order to gain a better 
understanding of the free parameters {ko,(T,p,n) that characterize the states. The first 
observable we shall consider is p^, that is a Dirac observable (and therefore a 'constant of 
the motion'). The operator is represented as p^ - x = ' Xi ci-nd its expectation value is 

thus given by 

(S><j>) 

Using Parseval's theorem we get 



Ixllphy 





dF 


J — oo 


dx 



dx',. 



2h 

I -^11 phy J — oo 



e\F{k)\^dk. 



(3.4) 



If we take our initial states F{k) = k^e ^of I ^^v^ then the integral takes the form 

2h 



{P'f>) 



Ixllphy 



(3.5) 



With the change of variables u = ^{k — fco), the expectation value (3.5) takes the form 



{P<f>) 



2h 



Ixllphy 



71. 



2n+3 ^oo 



U + 



V2ko 



a 



2n+2 



e"" du . 



(3.6) 



In the case of pure Gaussian states, namely for n = 0, we have 

,2 



{Pcf) 



hkn 



1 + 



(3.7) 



This is telling us that, for the approximation we are taking, namely fco ^ cr, then (p^) ^ hko. 
Therefore, as expected, the parameter ko is giving us a good measure of the expectation value 
of p^ which can be regarded as the conjugate variable to x, on the reduced phase space of 
the system (See Appendix |A]). Let us now compute the fluctuations of this operator. Let us 
start by computing 



iPl) 



2hH 



Ixllphy 



dF d^F 

dx dx^ 



-dx 



k''\F{k)\^dk', 



(3.8) 



phy -J —CO 
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For our generalized Gaussian states the integral takes the form 



(Pi) 



Which can be rewritten as, 



Ixllphy 



(3.9) 



Taking n = 0, 



(Pi) 



2hH 

l^llphy 



a 



V~2. 



iPl) 



2n+4 



u + 



a 



2n+3 



1 + 



h^kl 



Note that for large k^, that is for those states satisfying k^^ a, then (p^) 
Let us now compute the dispersion of the observable p^, {^P(j>Y = (P^) ~ (P^)^ 
n = case we have 



(3.10) 



For the 



{^P^f 



h^kl 



1 + 



3^2 
Akl 



h^kl 



1 + 



a 

4^J 



2^2 r 



1 - 



(3.11) 



which is a constant. Note that for k^ ^ a, the dispersion /S.p^ ^ ^. This is telling us that 
the parameter a has the interpretation one might have expected as the dispersion of the 
observable associated with the variable fc, which in this case corresponds to p^. 

Recall that the quantity that is conjugate to p<^, in the reduced model, corresponds 
to x(/3). Thus, it is natural to consider the family of operators x^, corresponding to the 
observable 'x at the time 0'. It is most natural to define the symmetric operator acting on 



initial states F, in the ^-representation the operator defined in Eq. (2.20) can be written as: 



xF 



id i 
dk^2k~ 



F. 



(3.12) 



The first two terms arise when Eq. (2.23) is replaced into Eq. (2.20). It is now straightfor- 



ward to compute the expectation value of the x operator on the generalized Gaussian states. 
It is easy to see that, for all values of n, 



(3.13) 



which confirms the expectation that the parameter —p represents the point in x space, where 
the Gaussian is peaked (when = 0). This also corresponds to the classical dynamics (as 
given by the effective theory) for the variable x (with our choice > 0). Furthermore, one 
can try to find the expectation value of the operator 



X 



id 
dk 



2k 



1 d 1 
kdk ^ 4P 



id i 
dk^2k 



+ 



Here it is important to note two important issues that arise when considering the operator 



defined by (3.12). The first one is that it involves a derivative operator and the second is the 
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factor in the operator. The derivative term imphes that one needs to consider carefully 
the boundary conditions at = 0, so that the boundary terms (when integrating by parts) 
does not contribute. In particular, we should have, for the x operator, for states defined by 
functions f{k) and g{k), that {g f k)\j^^Q = 0. This condition is, of course, satisfied by all 
of our states n > 0. The term 1/k imposes some fall-off conditions at the origin {k = 0) for 
the integrals to be finite. For instance, even when the expectation value of x is well defined 
in the n = Gaussians, these states are no longer in the domain of the operator. The state 
that results, when acting with the operator x on a n = state is not normalizable, and 
therefore, does not belong to the physical Hilbert space. We have to conclude that, in order 
to have the operator x in the set of observables we want to consider, we have to exclude the 
n = states and consider only n > 1. In the n = 1 case, it is straightforward (if lengthy) to 
compute the expectation value of the operator, 

= + ^ + 2P'#' + (3-14) 
" (l + tf) 

from which the fluctuations of the operator becomes, 

{Axf = {x% - {x)l = ) ^ (3.15) 

which also tells us that the parameter a has the expected interpretation of providing the 
inverse of the dispersion for the observable x. This is realized when ko ^ a and therefore 
the second term in the previous expression is very close to 1. 

We can now examine the uncertainty relations for the observables and x (for n = 1): 

4 ^+3-^' 



We note that, if ko ^ cr, then the uncertainty relations 

{Ap^){Ax) > ^ 

are very close to being saturated. As we shall see in the next section, these conditions are 
necessary for the generalized Gaussian states to represent acceptable semiclassical states. 



B. Volume Operator 



Let us start by considering the operator Vl^ig corresponding to the volume at 'time 



The expectation value of this operator is given by Eq. (2.15) so it suffices to compute the 
constants V± given by 



K 



dF 
dx 



dx - 



with K = An'jipY^ 
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which can be written as, 



V- 



_ 14 _ Ik 

^ Ti \\2 Ti \\2 

11X11 phy 11X11 phy ^ 



k{k T ia)F{k)F{k T ia) dk . 



(3.17) 



If we take the state as F{k) = /c"e ^'^ I ^i-pk ^ -^.j^g integral is 



± 



a 



This can be rewritten as, 



a 



\k{k T ia)]"+ie-('^-'=°) e-('=^'"-^«) dA; . (3.18) 
It is straightforward to see that, with the change of variable u = \/2k — \/2ko =i= the 



integral (3.18) can be written as 

n+l 



± 



V2a V4 
For the n = case we have, 



2u^ + AV2kou + Akl + 



1 n.+l 



e-'^^/'^^d^. 



± 



From which we get 



± 



4^2 



+ 4y2fcoM + Akl + e~"'/"'du , 



(3.19) 



If we take now a normalized states using Eq. (|3.3|), the result is 



V- 



± 



Ixllphy 



2a 



e e 



4A;o + 



+ 



(3.20) 



In order to calculate the dispersion of the volume operator, it is necessary to calculate the 
three integrals Wq, W+, W_. Let us start by computing the quantity Wq. 



W< 



ia2 







1 1 X 1 1 phy J —oo 

which can be written as 

Wn = 



dx"^ dx dx dx^ 



da;, with a2 



2^2 



Ixllphy 11x11 phy J — oo 



k^\F{k)\^dk. 



(3.21) 



(3.22) 



This is the same integral that the one in Eq. (3.8). Then, for a normalized Gaussian state 
with n = we only need to replace fi^ by a2 into Eq. (3.10) to obtain that 



Wo = 02^0 



1 + 



3a^" 
Akl_ 



(3.23) 
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Let us now compute the quantities W± = W±/||x||phy; with 



W± = ia^ / e 



dx"^ dx dx dx"^ 



3G 



dx, with as 



(3.24) 



that can be written as, 



W± = as k 



F{k){kT'2ia)F{kT'2ia) + {k±2ia)F{kT2ia)F{k) dk . (3.25) 
If we take the initial states F{k) = ^-i'^-^o)'^ / '^'^ ^'^vk ^ ^j^g integral is 

/CO 
-oo 

[(A; T 2ia)"+ie-(^^2ia-fco)V-^ + ± 2ia)"+ie-('=±2*a-fco)V-^]dA; (3.26) 
In order to do the integral 

j'OO 

/(±) = / p+2e-('=-'=o)V-^ ± 2m)"+^e-('=^2ia-fco)V-2d^ _ 



It is straightforward to show that with the change of variables u = \/2{k — k^ ±ia), we get 

n+2 / \ n+1 



U 

71 



+ fco =F 



u 
71 



72 



.Jla^la"^ /-oo 



(^ + ^o) (^^ + A;oT^a)e-"'/'^'dM. (3.27) 



The real part of /(±) is 

p2c?ja'^ poo 

Me[/(±)] : 



n+1 



71 



+ A;o ) e-"'/"'dM. (3.28) 



Then Me[/(+)] = ]Re[/(— )], and the imaginary part of /(±) is 



Im[/(±)] = Ta- 



{- 



^kA 



n+1 



theref ore, _ Im [/(+)] = — Im[/(— )]. Using the real and imaginary parts, we have that the 
integral W± takes the form 

W± = a3e^'^"[/(+) + /(-)] = 2«3e^'^"Me[J(+)] . (3.29) 

This tells us that W± is real valued (as expected) and that it is only necessary to calculate 
the integral (3.28). Taking n = the integral (3.28) takes the form 

2 



ie[I{±)] 



,2aV<T^ 



U 



x/2 



V~2 



u 



V2 



+ fco e-" dw 



22/3 



3y A;o + ^(3A;2 + a^) + (A;^ + a^k^^ e'^^'^^'du . 
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Using the integrals from Appendix [B| the integral is 



ie[I{±)] = -koa^ + kla + a^koa 



V2 



from which we have that, 



Including now the normalization of the state we have 



kl + ^(T^ + 
u 4 



(3.30) 



(3.31) 



(3.32) 



From the last equation and Eq. (3.20) we can observe that the difference in the asymptotic 



relative fluctuations after and before the bounce, as defined by Eq. (2.18), vanishes. That 



is, D = 0. It is important to note that this is a general result for any Gaussian state of the 



form of Eq. (3.1), for all values of n = 0, 1, 2, as was expected to happen [T8] . 



IV. SEMICLASSICALITY CONDITIONS 



So far we have considered generalized Gaussian states as initial states for physical states 
in the exactly solvable fc=0 LQC model with a massless scalar field. We have seen that 
the parameters that define the states, in the pure Gaussian case with n = 0, have the 
expected interpretation: The parameter k^ is related to the expectation value of p<^, —p to 
the expectation value of £, and ha /2 to the dispersion of p^. The only condition that we 
have imposed so far is k^ S> cr, which guaranties the validity of our approximation. We 
have also seen that, for states which satisfy this consistency condition, there are at least two 
results that follow. First, it warranties that the relative dispersion of p^ is small. Second, 
this conditions implies that the uncertainty relations become saturated in the sense that 
(Ap^)(Ax) ^n/2. 

From this perspective, it could seem that any value of k^ and a, provided they satisfy 
fco ^ cr, might be acceptable to define semiclassical states. The purpose of this section is 
to explore this issue further and answer the following question: Are there more stringent 
conditions that one must impose in order to have semiclassical states that will further restrict 
the possible values of /cq and cr? As we shall see, the answer is in the affirmative. 



A. Asymptotic Volume 

In order to answer this question we shall first consider the volume operator. Since all 
physical states have the property that the expectation value of the volume operator V^|<^ at 
time follow the same functional form and, therefore, follow for large the same dynamics of 
the classical dynamics (i.e. Einstein's equations) ^(0) ~ e^""^, we need more criteria to select 
those states that are semiclassical. The obvious strategy is to consider the state's relative 

dispersion ^(AV^)/(y)^ of the volume operator at time (p. One expects that semiclassical 

states will have a very small relative dispersion 'at late times' when the dynamics approaches 
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the classical dynamics. It is then natural to consider the asymptotic relative dispersions given 

by 

A,:= lim f^V = H^-l. (4.1) 

It is straightforward to find the analytical expression for these quantities in our states (in 
the case n = 0) to be 



A± = - 1 = e°^/-^ \ :V - 1 . (4.2) 



It is interesting to note that, in the previous expression, there is some competition between 

2 / 2 m 

the factor e" ''^ and the second part. The exponential becomes close to one if a is large 
compared to a, but can become very large if a becomes too small. On the other hand, the 
larger a becomes, the larger the second term. Thus, that there must be an optimal value 
a for which A-t is the smallest. The form of A± can be seen in Fig. [TJ Another important 
aspect to consider is the following. So far, we have not imposed any condition on p^. That 
is, we have not said that it has to be 'large', mainly because in the classical theory there 
is no dimension- full quantity with respect to which one can compare it. In the quantum 



theory, with the introduction of h we do have such scale. As we can see from Eq.(4.2) and 



the expression for (p,/,), the quantity ha has the same dimensions of p^. Furthermore, if we 



now impose the condition ^ a, all the terms in parenthesis in Eq.(4.2) are small, so we 
can indeed try to find the optimum value for a. Note also that the condition a implies 
{Vtt)) ^ ^ct- It is in this sense that we can call the momentum 'large'. 



The condition that the expression (4.2) be an extrema is a cubical equation, and the only 



physically interesting solution can be approximated as, 

a'^ ^2akQ + '^a^ . (4.3) 

Since we require that k^^ a, then the second term is very small compared to the first one 
and we can approximate a ~ a/2 a k^. It is important to stress that this condition selects 
a preferred value for a for which the state is most semiclassical. As can be seen from the 
Figure [T| if we make a slightly smaller (in an attempt to 'make the dispersion small') the 
asymptotic relative dispersion in volume becomes very large, making the state not a good 
candidate for a semiclassical state. 

If we now introduce this value for a into the relative fiuctuation function, then this can 
be approximated by 



lim 

^—^■±00 



(AlO 

{V? 



(.0) - I (4.4) 



This approximation has an error of 0.0002% when k^ = 30000 (in Planck units), and becomes 
better as ko grows. This tells us the precise relation between the parameter ko that controls 
the expectation value of p,/, and the minimum value of the asymptotic relative fiuctuation of 
the volume operator. The a dependence of the asymptotic relative dispersion of the volume 
can be seen in Figure [TJ with k^ = 30000. Two features characterize this plot, namely, 
the first one is that for < the function is exponential so it grows very quickly as 
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FIG. 1: The asymptotic relative fluctuation A-t is plotted as a function of with fco = 30000 in 
Planck units, where a ~ 6.14. Here we have chosen n = and p = 0. Two features characterize 
the behavior of the function near the minimum a. The first one is that for < a"^ the function is 
exponential; the second one is that for o"^ > o"^ the function has a polynomial behavior. 

one approaches zero; the second feature is that for cr^ > the function has a polynomial 
behavior, so one can have values of a > o" without A± changing too much. 

To summarize, we have seen that not any Gaussian state will be an admissible semi- 
classical state. In order to minimize the asymptotic relative fluctuation A±, one has to 
choose the parameter a carefully as a function of the other relevant parameter (in this case, 
the parameter ko). That is, if we choose a to be of the order of o"^ ~ 2akQ, we can then 
minimize the asymptotic relative dispersion of volume after and before the bounce (which 
are the same, as we have seen). Gaussian states for which this condition is not satisfied will 
have, in general, large quantum fluctuations in the regime where we expect them to be small 
-for large times after and before the bounce- so they can not be regarded as semiclassical. 

Let us then assume that we have chosen the value a for a. The question we want to 
ask now is the following: Can we reproduce some of the features that are characteristic of 
the so called "effective dynamics"? As we have discussed before, the 'classical' Hamiltonian 
constraint one obtains from replacing classical quantities like connection and curvature by 
a corresponding holonomy functions with a finite parameter A, can be seen as either a 
starting point for 'polymer quantization', or as the 'classical limit' of a loop quantized 
theory. In any case, we expect that the semiclassical states of the quantum theory approach 
this 'classical theory' in the appropriate regime. In other words, we want to know when 
the effective dynamics is a good approximation to the full quantum dynamics. Is there a 
choice of parameters in the Gaussian states for which the effective dynamics is not a good 
approximation? We shall answer these questions in the remainder of this section. 

B. How Semiclassical is the Bounce? 

Another important issue is the question of how semiclassical the state is at the bounce. 
For instance, if the state has a very small relative asymptotic dispersion A± at late times, 
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one might want to know what the relative dispersion is at the bounce and how it compares to 
A-I-. A plot of the relative fluctuation, as a function of internal time can be seen in Fig. |2| 
where it is shown that the relative fluctuation is bounded, it attains its minimum value 
at the bounce and approaches A± asymptotically. Thus, somewhat surprisingly, the state 
does not become very "quantum-like" at the Planck scale but rather preserves its coherence 
across the bounce. 

Let us now try to estimate the value of the relative dispersion in comparison to its 



asymptotic value. We know that, at the bounce, 4 
2y/V^, and 

iAV)l^ =Wo- 4Wr_ 
Then the relative fluctuation in the bounce is 



fin 

2a 



Aa 



In 



W- 



w+ 



then {V)^ 



Wo + 2.JW-W^ 



4-lA VI 



1. 



(4.5) 



(4.6) 



For n = in the Gaussian states we have 



1 {AVf 



2 lyy 



<f>^±oo 2(1 + 



(4.7) 



Now, if we impose fco ^ o" ^ a, which is satisfled if we choose a = a = y/2ako, then the 




FIG. 2: The relative fluctuation {AVf/{Vf is plotted as a function of internal time (p with 
ko = 30000, in Planck units where a ~ 6.14. Note that the minimum of this quantity is at the 
bounce (here corresponding to i;^> = 0). It is found that the value at the bounce is approximately 
1/2 of its asymptotic value for the generalized Gaussian states. 



last equation has the form. 



1 {AVf 
2 



0{a'/a^) + 0{<yykl) 



¥±00 
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This tells us that when ^ a ^ a, the relative volume fluctuation at the bounce becomes 



{Avy 



(AVf 



2 ^ 



(4.8) 



>itoo 



This approximation, for instance, has an error of 0.007% when ko = 30000, and it becomes 
smaller as ko grows. If we use the Eqs. (4.8) and (4.4) and choose the value of a that 
minimized the relative fluctuation we get 



(AV^)^ 



a 

2ko 



This approximation have an error of 0.01% when ko = 30000, and becomes smaller as we 
increase ko. We have performed intensive numerical explorations and have seen that, for 
n > 1, the relative fluctuation is again smaller at the bounce and approaches 1/2 of its 
asymptotic value. 

Thus, we can conclude that states that satisfy the conditions fco ^ ^ a, imposed by 
our requirements on the asymptotic relative dispersion of the volume through the choice 
a = a, exhibit also semiclassical behavior at the bounce, in the sense that the relative 
dispersion in volume is bounded and small. Even more, the relative fluctuation at the 
bounce is approximately 1/2 of the asymptotic value. These analytical results give support 
to the numerical explorations reported in [5l |6] . 



C. Volume and Energy Density at the Bounce 

Having shown that a semiclassical quantum state behaves also semiclassically at the 
bounce, it is then natural to ask whether the state yields observables -through their ex- 
pectation values- that follow some classical trajectories. It turns out that it does, but the 
'classical theory' it approaches is not the standard classical theory (GR+ massless scalar 
field), but rather an effective theory, as defined by an effective Hamiltonian constraint C\ 
|23j . which contains the quantum geometry scale A. The effective dynamics, is generated by 
the effective Hamiltonian constraint 

and has several important features. The first one is that all trajectories have a bounce that 
occurs at the same critical density pcrit = snc-y^x'^ • "^^^ volume reaches a minimum value at 
the bounce that depends on the observable p<^: 



K.in = y ^P<^- (4.9) 

In the quantum theory, we have seen that every physical states has a bounce where the 
minimum volume is given by Vbounce = 2a/V^ V_. It is then natural to find the corresponding 
volume at the bounce for our family of Gaussian states. It is straightforward to see that the 
minimal value of the volume, for the case n = 0, is 



T/ _ ' "-Pi'' pOi- rid- 7 

•^bounce ^ n,Q 

a 



1 + 



(4.10) 
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which is of the form, 



^bounce = ko + 0{a'/a') + Oia'/kl) + 0{a^/kl] 

a 

If we recall our previous result that (p,^) = hko + 0{a'^/k'^), and assume the semiclassically 
conditions found in previous sections, namely that a ^ a and /cq ^ o" (as is the case when 



a = a =~ \^2ako and ko ^ a) in Eq. (4.10), we can compare the last two equations to see 



bounce 



Using the value a = ^/ITttG, we conclude that 



bounce 



Klin, (4.11) 



which is what we wanted to show. Note that the higher the value of ko, the better the 
approximation becomes. In other words, if ko were not large enough (compared to a), not 
only would the asymptotic relative fluctuation of volume be large, but the effective equations 
would fail to be a good approximation to the quantum dynamics. This seems to indicate that 
quantum Gaussian states with a higher value of ko behave more classically, in contradiction 
with the classical intuition that tells us that 'rescaling of fco' is physically irrelevant^. Thus, 
we have to conclude that this rescaling symmetry of the classical theory is broken in the 
quantum theory. Note however, that even when the exact symmetry is broken, one might 
expect to regain it in the limit of large ko which, as we have seen, corresponds to (p^) ^ ha. 
For an in-depth discussion regarding this scaling freedom see |26] . 

Another important result for this solvable model is that the energy density p is absolutely 
bounded by the critical density pcrit ilZj- It is then natural to ask what the behavior of 
the energy density at the bounce is, for semi classical states. One might imagine that, for 
instance, the more semiclassical the state, the higher the density at the bounce. In this part 
we shall explore this question by taking several quantities to measure density. For instance. 



Ml 

2{V)- 

straightforward to find this quantity for the n = states. 



the simplest one would be p = a quantity shown to be bounded by pcrit It is 



e-V-^ 1 + 4 

P = Pent — ^ ^ (4.12) 

Note that p < pcrit and one approaches it as /cq — oo, assuming again that we have chosen 
cr ~ yj2a ko and ko ^ a. One can also see that, if we were to fix ko and leave a free, the 
density grows and approaches pcnt as a grows. Thus, asking for the density to approach 
the critical density is not a very stringent condition on a as the relative fluctuations in 
volume was. Note also that, if ko was small enough to be of the order of a, the density at 
the bounce would be far from the critical density giving yet another indication that those 
states are not semiclassical. The other quantities representing density that one can build, by 
taking for instance (p^) and (t^^), give expressions that have the same qualitative behavior, 
approaching pcrit as ko grows. 



^ Recall that the k—Q FRW model has a rescaling symmetry of the equations of motion, and of the underlying 
spacetime metric given by {V, /?, (j),Prt,) — >■ {(-V, (3, cf), £p^) for £ a constant. This can be understood as coming 
from the freedom of choosing arbitrary fiducial cells for the formulation of the theory [22] • 
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D. Where is the Classical Region? 

0.75 r 



0.70 



0.65 
0.60 
0.55 



0.50 



0.45 



100000 200000 300000 400000 500000 

FIG. 3: The values 4> (p for which one reaches the classical region are plotted as a function of 
ko, where n=0, a = a and p = at (j) = 0. In the plot, 6 = 1 dotted, 5 = 0.5 dashed, S = 0.25 
dotdashed, S = 0.125 down line, 6 = 0.0625 middle line, and 6 = 0.03125 the top line. Note that 
for high values of fco, the time (p remains a constant. 



p 

Peril 

0.020 - - 
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100000 200000 300000 400000 500000 

FIG. 4: Values of density are plotted for different values of 5((/>), as function of ko, where n=0, 
a = a and p = at (j) = 0. Here, /01234/Pcrit for which 5 = 1 dotted, S = 0.5 dashed, 5 = 0.25 
dotdashed, 6 = 0.125 up line, 6 = 0.0625 middle line, and 6 = 0.03125 the bottom line. 

In the LQC description of the dynamics as governed by the effective Hamiltonian, one 
can unambiguously ask when one recovers the GR classical dynamics. One way of measuring 
this could be in terms of the variable P (when A/3 <^ 1), or in terms of the energy density 
(p ^ Pcrit)- However if one wants to ask a similar question in the quantum realm, one 
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immediately faces the following dilemma. As we have seen in the previous part, a Gaussian 
state remains semiclassical across and after the bounce. In fact, the closer to the bounce, the 
more semiclassical it is. We can not therefore ask when does the state become semiclassical, 
based on the behavior of fluctuations of the relevant operators. Still, we would like to have 
the notion of a transition from a quantum gravity dominated regime to a the regime where 
the quantum gravity effects are small and one is in the "classical regime" . 

For this purpose, we want to propose such a criteria, and define the (percentual) error 
between the relative fluctuation and the asymptotic relative fluctuation, 



5{<p) := ^ X 100 . (4.13) 

We expect that this error will be small when we are approaching the classical region. Let us 
now motivate this proposal. While we can not define a transition in terms of a change in the 
relative fiuctuations of the basic operators, we could compare them to those of the standard 
'Wheeler DeWitt' model. For we know that the expectation values in the WDW theory 
follow closely the classical dynamics of general relativity. There is indeed a canonical way 
of mapping a LQC state to a WDW state, so the comparison of expectation values is well 
defined |l7j. We have done that for our Gaussian states and computed the corresponding 
relative fiuctuation for the volume. It turns out that, for the WDW dynamics, the relative 
fiuctuation is constant in evolution, and corresponds precisely to the asymptotic relative 
dispersion of LQC. Thus, the LQC evolution not only approximates the WDW one in terms 
of expectation values but also in terms of the relative fiuctuations. It is then justified to 
measure the transition to the 'classical regime' in terms of how close the LQC evolution is 
to the corresponding WDW dynamics. 

In order to study where the classical region begins, we found the value of (j) for which 
(5 = 1/2™" with m = 0, 5, taking = at the bounce. The values of that satisfy these 

conditions are plotted in figure jsj Then we evaluated the densities pi = , p2 = 

' Ps ~ oTife ' P4 ~ ^uv2\ ' this values of 0. The density is plotted in figure 4 This plot 
does not depend of the used density {pi, P2, Ps, Pa)- From figured we can see that 



^M.M. (4.14) 

Pcrit 5 



for /co ^ a, cr = a = y/2akQ, p = and n = 0. These results are similar if n > 0. 
Furthermore, if we choose p ^ 0, the only change is in the value of at the bounce. If we 
change the value of a for a fixed ko, then for larger (than a) values, the time of arrival to 
the classical region is smaller. In that case, the value of the density is closer to the critical 
one, and in the limit of large a the classical region is 'at the bounce', which does not make 
much sense. In this way the need to bound a to be near the value a is manifested. 



The relation given by Eq. (4.14) is telling us the order of the quantum corrections. Thus, 
if we say that the quantum corrections are negligible when, say, 6 ~ 10~^° then the density 
will be of the order p/pcrit ~ 2 x 10~^^. Recall that Gaussian states are symmetric about 
the bounce, namely A+ = A_, therefore all the conclusions of this part equally apply before 
the bounce. 

Let us now summarize this section. We have seen that appropriate physical conditions 
imposed on the behavior of basic physical observables are enough to constraint the param- 



eters that characterize the generalized Gaussian states (3.1), and obtain truly semiclassical 
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states. Of the three continuous parameters {ko, a,p) that characterize the states, two of them 
are fixed by semiclassicahty considerations, while the choice of p represents a true choice of 
initial condition. We have seen that asking that the relative fluctuations of the observable 

be small imposes that ^ a. Smallness of the asymptotic relative fluctuation of the 
volume implies then that kQ ^ a ^ a. If these conditions are satisfied, we saw that the 
semiclassical states at 'late times' remain semiclassical across the bounce and that the vol- 
ume and density at the bounce are very well approximated by the effective theory. We have 
also seen that the scaling symmetry of the classical theory is not present, within the class 
of states under consideration, when ko is of the order of a. However, in the limit of large 
kf), all the properties of the state remain invariant -in terms of being well approximated 
by the classical effective theory-. This leads us to conclude that the scaling symmetry is 
approximately recovered for "large p^" (for details regarding this issue, see [26j). 

Still, the Gaussian states we have considered are not the most general states one can 
consider. In the following section we shall explore the so called squeezed states, a general- 
ization of the Gaussian states and compare their semiclassicahty properties to those of the 
Gaussian states. 



V. SQUEEZED STATES 

In this section we consider states that generalize the Gaussian states considered so far. 
The strategy will be to consider squeezed states that are nearby the Gaussian semiclassical 
states and compare their properties. In particular, we would like to know if the squeezed 
states can improve the samiclassicality properties such a better behavior of relative fluctu- 
ations of volume or a better approximation to the 'effective theory'. 

The generalized Gaussian state we have considered so-far had three free (real) parameters 
(p, ko, a) and one discrete parameter (n). We saw that, given a point of the physical phase 
space, we could approximate the internal dynamics of the theory by a choice of p, and that 
semiclassicahty imposes conditions on the two other parameters ko and a. 

Let us now take the initial states as 

f j^n^-nik-p)^ forA;>0, r?, /3 G C 
F(k) = l ' (5.1) 

[ 0, for A; < 0, 

depending on two complex parameters (?7,/3). One can reduce to a Gaussian state from a 
squeezed state by setting 

Im(r/) =: r]i = 0, Me(r/) =: r]R = ^, Me(/3) =: = ko, 2r]Rf3i = p . (5.2) 

Thus, we see that the extra parameter in the definition of the states if given by the imaginary 
part of 7]. The norm of this state in k space is 



oo 



(X,x)phy = 2 / ke-e-^^'-^^\-^^'-f^^'dk. 

J — oo 

Using relation (B9) for t = the last integral can be written as 

/oo 
-oo 
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where 



a = I3r- I3j 



Vr 



(5.4) 



If we want the convergence of the integral then rjj^ > and if we want to select the positive 
frequency then a ^ i.e. (3r — f3i^ ^ With this assumptions the integral is well 
defined, and the error is small, as shown in Appendix [Cj 

If we take n = and use the change of variables u = k — a we have 



Ixllphy 



/oo 
{u + a)e-'^^''"dk. 
-oo 



Then 



Ixllphy 




(5.5) 



(5.6) 



A. Elementary Observables 

We can compute the expectation value and dispersion of the fundamental observable 



2n 



Ixllphy 



If we use Eq. (B9), with t = 0, then we obtain 

2h 



If we take = 0, 



Ixllphy 



2h 



I 1 1 phy 



TC 



2vr 



This becomes, using Eq. (5.6) 



h 

a 



Now we want to calculate the expectation value of p| 

2^2 rco 



k'\F{k)\'dk. 



phy — oo 



In the squeezed states (5.1) we get 

iPl) 



I ^11 phy 



pn+3g-r)(fc-/3)2 g-r,(fc-/3)2 



If we now use Eq. (5.6), the result is for n = 

iPl) = 



+ a" 



(5.7) 



(5.8) 



(5.9) 



(5.10) 



(5.11) 



(5.12) 



(5.13) 
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From this we can find the dispersion and get, 



Ar]R 16r7|a2 
Writing this exphcitly in the original variables is 



(5.14) 



(5.15) 



B. Volume 

We shall compute the quantities V± and Wo,± from which one can find the expectation 
value and dispersion of the volume operator. The relevant coefficients are 



K 

a 

K 

a 



i[k ^ ia] [k T mj^e-^^^'^'^-^^'dA; 



■oo 
oo 



(5.16) 



with K = AnjipiX. As in previous subsections we will only write down explicitly the n = 
case. This expression can be put on the form, 



V± = — exp 
aa 



a 



1 bl 

h — H 

8vr 2 8 



(5.17) 



where fe± = a=Ff^ = /3ij— ^ ± f)- In the Gaussian variables it takes the form 



V- 



± 



- ko-7]j — 

a \ 2 



exp 



(^±pa)(l + ,„V) 



X 



1/ /2 \\2 2 

a 1 / , 2 1 P'^ I a \ \ ct 



± 

2 2 



+ 



(5.18) 



which reduces to the Gaussian case when rji = 0. 

The next quantities we need to compute in order to find the dispersion of the volume 
operator are the coefficients (Wq, W±). The first one is given by 



Wo = 2a2 / k^\F{k)\'^dk, with as 



3G 



(5.19) 



This is the same integral that the one in Eq. (5.11). Then, in order to consider the proper 
normalization for the squeezed state with n = 0, we only need to replace ti^ by 0^2 into Eq. 



(5.13), the result is 



Wq = a2 



^Vr V Vr 



(5.20) 
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The other quantities are given by the expression, 



W± = as k 



F{k){k T 2ia)F{k T 2m) + C.C. 



dk. 



2a3Me [I±] 



(5.21) 



with 03 



and I± = f^^k'^F{k){k =1= 2ia)F{k =1= 2ia)dk. This integral can be 
performed for n = and we get 



W± = as exp 



(-1)1 


X 







1 



2?7^a 4r]ji 



+ b2 + o? 



(5.22) 



where = a^a^ = Pr — ^ (/?/ ±0;). It is now straightforward to find A-t, the asymptotic 
relative dispersion for volume, 



A 



± 



exp 



2 

Tr 



^ 1 



1 I fc± I a2 
8r;ji 2 8 



-1', 



(5.23) 



Let us now explore those conditions that select semiclassical states. 



C. Semiclassicality Conditions 

In this part, we want to see how the semiclassicality conditions are satisfied if we move 
along the 'squeezing parameter' rjf, off the Gaussian semiclassical states, corresponding to 
rji = 0. Instead of performing an exhaustive analysis of this question, as we have done in 
previous sections for Gaussian states, we have computed the relevant quantities and plotted 
them. 

Recall that for Gaussian states, asking that the asymptotic relative dispersion in volume 
be a minimum, selected an optimal 'width' a = a ~ \/2ak^ for the Gaussian. This value 
had also the property that the density at the bounce was very close to the critical density 
Pcrit- In Fig. [5] we have plotted the asymptotic relative dispersion and density at the bounce 
as functions of a and 777. The rjj = slice corresponds then to the values found for Gaussian 
states. It can not be appreciated in the figure but, for a fixed a, the asymptotic relative 
dispersion A±, as a function of rji does not attain its minimum at zero, but rather at some 
value fjj ^ for A+ (and rjj = —fji for A_). Thus, if one chooses, for the parameter of 
the squeezed state the value r]i = fji, one is 'improving' the asymptotic relative dispersion 
after the bounce, and, at the same time, increasing the value of the asymptotic dispersion 
before the bounce. In this process, one is introducing an asymmetry in the fluctuations. 
This fact has been used to suggest the possibility, for example, that the asymmetry could 
be so large to spoil the semiclassicality properties across the bounce. As we shall see below, 
this possibility is however, not realized. What one can indeed see from the figure is that, as 
we move further away from the Gaussian states, that is, as we increase to values |?7/| > {fjil, 
the asymptotic relative dispersion increases steeply, for both signs of 77/. As can be seen 
from the figure, the density at the bounce (in blue) decreases as we go away from rji = in 
both directions, reaching densities much less than the critical value very fast. 
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FIG. 5: Density and asymptotic relative fluctuation A+. In the plot = o"^ : (0,100000) , 

r// : (-0.001,0.001) , /Sr = ko = 30000 , /3/ = and n = 0. The maximal value of the density 
correspond to the Gaussian states. Is clear that the squeezed states have a maximal density close 
to pcrit only if the state is near to the Gaussian states. 




FIG. 6: Volume at the Bounce. In the plot l/rjR = : {0,8 x 10^) , r]j : (-0.0001,0.0001) , 
I^R = ho = 30000 , /?/ = and n = 0. The global minimum is over the Gaussian states and 
has an exponential behavior if we move far away from the Gaussian state or if cj^ is smaller that 
= 2ako. 
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Both phenomena suggest that we can indeed consider squeezed states as semiclassical 
states provided we remain very close to the Gaussian states. The question is then how close 
we can be in order to still have semiclassical behavior as exhibited by the Gaussian states. 
For this purpose, we have plotted in Fig. [6] the volume at the bounce as function of a and 
r]i. We again see that, as was the case for the density at the bounce, the volume at the 
bounce has indeed a minimum at the Gaussian states, namely, for rjj = 0. Using both 
relative volume and density at the bounce, we can then define an interval, as function of a, 
in which the parameter rjj can take values and the state still be considered as semiclassical. 
More precisely, we define some tolerance for the asymptotic relative dispersion A+ = A, and 
find the maximum value of r^™^^ for which the relative dispersion is below the value A. We 
then find the value of the density at the bounce for the value 7]^^^- We have made extensive 
explorations for different values of A and ko, and have found that, for a fixed value of A, 
the larger the value of ko, the smaller the allowed interval in 77/ . For example, for A = 0.01, 
the dependence of rjf^^ on ko can be approximated (for large ko, a = a, n = and p = 0) 
as r^f ~ l /v/47000A ;o. The value of pb at that point is equal to 0.99pcrit- (For A = 0.001, 
rif^'' ~ l/^/mmO~k) and pb = 0.999pcrif) 




FIG. 7: Relative error of the difference in the asymptotic values of the relative fluctuation . 

In the plot 1/rjR = cj^ : (0, 10^) , r// : (-0.0003,0.0003) , /3r = feo = 30000 , ^7 = and n = 0. 
Outside of the vertical one the squeezed state have an asymptotic relative fluctuations (A_|_ and 
A_) bigger that 1. One line mark the Gaussian states {r]j = 0) that is just when the curve is zero 
and the other line is for states with o"^ = 2ako. 

Let us now return to the issue of the asymmetry across the bounce of the relative disper- 
sion in volume. As we have mentioned, without any control on the value of the parameters, 
one could in principle have a large asymmetry that translates into a loss of semiclassicality 
across the bounce. This scenario was suggested in [21] , but where the analysis was made for 
an effective formalism not directly derived from LQC and with little control on the semi- 
classicality conditions. Here we are better equipped to make precise statements about this 
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asymmetry based on the control we have on analytical expressions and the possible values 
the parameters of the states can have. In Fig. [7] we have plotted the relative difference 
A+/A_ — 1 in asymptotic relative dispersion, as a function of a and r/j. We have included, 
as vertical surfaces, the boundary of the allowed region in parameter space. Two slices in 
the figure are worth mentioning. The curve ?7/ = 0, which is viewed almost vertical in the 
figure corresponds to the Gaussian states. On that curve, A+ = A_ so the relative differ- 
ence vanishes. The second highlighted curve corresponds to the choice (7 = 5", which is the 
extension to r// 7^ of the optimal Gaussian states. Note that this curve starts at zero, 
increases (for positive r//), reaches a maximum and then decreases tending to zero. The 
curve is antisymmetric with respect to r//. Thus, there are two points for which |A_|_/A_ — 1| 
reaches a maximum. What we have seen, by taking many different values of fco? is that this 
maximum value is indeed, extremely small. For instance, for = 10^, the maximum is 
about 10~^. As we increase /cq this value decreases even further. These results show that 
indeed, semiclassical states are very symmetric and that semiclassicality is preserved across 
the bounce. One should also note that this quantitative analysis invalidates the claims of 
[20] . and [21] regarding the allowed asymmetry in fluctuations across the bounce and 'cosmic 
forgetfulness', and supports the results of [18] and [IH]^. For instance, in the first reference 
of [21] , it is claimed that the relative change in relative asymptotic dispersion in volume can 
be as large as 20 for semiclassical states, while we have demonstrated that this quantity is 
several orders of magnitude smaller for realistic values of the parameter (of the order of 
10~^ for ko = 10^, and much smaller for the values expected to represent realistic universes 
{ko ~ 10^°*^ or larger)). 

VI. DISCUSSION 

Let us summarize our results. We have defined coherent Gaussian states as candidates 
for semiclassical states 'peaked' around points of the classical (physical) phase space. We 
have imposed consistency conditions on such states by asking that the relative dispersion 
of the volume be small at late times after the bounce. This condition implies that one of 
the parameters, yielding the momentum of the scalar field, be large and fixes one of the 
other parameters. We can therefore find a canonical Gaussian state for a given point of the 
classical phase space corresponding to a classical solution at late times. When exploring 
the properties of these states near the quantum bounce we found that they behave also 
semiclassically in the deep quantum region; the relative fluctuations of volume are of the 
same order and smaller than in the asymptotic, large volume, region. Furthermore, for those 
states that have small fluctuations, the expectation value of the volume and energy density at 
the bounce are very well approximated by the so called effective theory, thus giving support 
to the claim that this theory is a good approximation to the quantum dynamics, even in the 
deep Planck regime. Next, we introduced squeezed states, a one parameter generalization 
of Gaussian states and studied their properties in a vicinity of the semiclassical Gaussian 



One should recall that the model defined in [20j . even though using a quantization quite different than 
that in [6j could be seen, in a loose sense, as a different factor ordering from the choice made in LQC and 
here. For a discussion on the differences between the assumptions made in the model of [20, and standard 
LQC see [T^, and for a detailed criticism of some claims made in [20], see [TH]- Our results here can be 
seen as giving further validity to the arguments presented in flB] against the claims of [Ml [2T] . 
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states. As we showed, the range of this extra parameter is severely restricted if we want to 
maintain a small relative fluctuation for the volume. Furthermore, if one departs from the 
Gaussian states too much, both the volume and density at the bounce differ very rapidly 
from the value on the Gaussian states (and also the 'effective theory'). Thus, we are lead 
to conclude that the Gaussian states exhibit a better semiclassical behavior than squeezed 
states. 

An important issue that can be studied quantitatively in this solvable model is that of the 
asymmetry on the volume fluctuations. On Gaussian states the fluctuations are the same 
before and after the bounce. As one introduces the squeezing parameter, the fluctuations 
are no longer symmetric. In fact, depending on the choice of sign of the parameter, either 
one of the fluctuations becomes smaller and reaches an absolute minimum in the vicinity of 
the Gaussian states. At the same time, the relative fluctuation on the opposite side of the 
bounce increases. The question is then how big can this asymmetry be? What we saw is 
that for a given, somewhat arbitrary choice of parameters, that is still far from those giving 
rise to a realistic 'large' universe, the relative difference in relative dispersion is very small 
(of the order of 10~^), and decreases for more realistic values of the parameters. Thus, even 
if we were to choose the squeezing parameter such that this asymmetry is maximized, the 
relative change is so small that, for all practical purposes, the state is symmetric. A state 
that is semiclassical on one side of the bounce not only remains semiclassical on the other 
side, but maintains its coherence. This quantitative analysis therefore invalidates claims of 
'loss of coherence' across the bounce where it was claimed that the relative change in 
(relative) dispersion can be several orders of magnitude larger. 

One particular feature of the classical description of the system under consideration is 
that it possesses a rescaling symmetry. In the spacetime description, we can rescale the 
volume by a constant and the physical properties of the spacetime remain invariant. In 
the phase space description this is manifested by a symmetry in the equations of motion. 
This symmetry, however, can not be treated as gauge, in the same sense that Hamiltonian 
symmetries are. One can understand the origin of this symmetry by recalling that the k=0 
model needs a fiducial volume for its Hamiltonian description, and this choice is completely 
arbitrary. In turn, this means that there is no scale with respect to which one could compare 
the momentum of the scalar field; there is no meaning to the statement that ^p^ is large'. 
As we have seen, in the quantum theory Planck constant introduces such a scale and this is 
manifested in the behavior of certain quantum observables. One of the semiclassicality con- 
ditions we have found is that the larger the expectation value of (p^), the more semiclassical 
the state is. What we have also seen is that the rescaling symmetry is only approximately 
recovered, for our class of states, in the 'very large (p^) limit', which can be taken as the 
classical limit of the system. For further details on this issue, see [26] . 
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Appendix A: Effective Reduced Phase Space 



The kinematical phase space of the model is F = with coordinates V, (f),p^)- Here 
(f) is the matter scalar field, its canonical momentum, V the volume (of the fiducial cell) 
and P its canonical momentum. The effective Hamiltonian constraint of the system is [23] : 

which is labeled by a real parameter A with dimensions of length. The symplectic structure 
is given by: 

^^0 = T^dV A d/3 + dp^ A d0 (Al) 
47rG7 

which induces the following Poisson brackets 

{/3,y}o = 47rG7 , {0,Mo = l- 

The Hamiltonian constraint defines a hypersurface T in phase space. It can be taken as a 
disjoint union of hypersuperfaces where /3 G (n|, {n + 1)^) with n G Z. In what follows we 
shall only consider the connected component where /3 G (—^,0). Note that the points where 
/9 = n|, n G Z are not in the constrained surface. These are degenerate points (hyperplanes) 
where the system does not have a well defined dynamics. We shall see that physically, these 
points should be excluded. 

We shall now consider the pullback Cl of the symplectic structure to the hypersurface T. 
That is, we find Cl such that Qq = t * Q, for t : F — )■ F the embedding of F into F. 

Let us write 



then, on the constrained surface we have 

p^ ^ ±r]V sin(A/3) (A2) 

The pullback of the gradient of Hx to F yields, 

r]'^2VdV sin(A/3)2 + 27]^V^ sin(A/3) cos(A/3)Ad/3 ^ 2p^dp^ (A3) 
±r7[sin(A6)d\/ + V cos(A/3) Ad/3] ^ dp^ (A4) 

which can be used for solving dp<^, dV or d/3 and inserting back in Eq. (Al). Each of this 
choices will depend on our choice of coordinates for the constrained surface. Let us study 
two of such parametrizations. In the first one we solve for p^, which means that we use 
{V,(3,(l)) as coordinates. Then 

fii = -j^dV A d/3 ± r][sm{Xf3)dV + V cos(A/3)Ad/3] A d0 
47rG7 

which gives us the expression for the two disconnected hypersuperfaces parametrized by 
positive and negative values for p^f,. The value p,/, = is excluded if (3 rij orV^T^O, with 
n G Z. If we take > we get 

Qi = -^dV Ad^ + T] sin(A/3)d1/ Ad4> + 7]XV cos(A/3) d/5 A # 
47rG7 
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If we now solve for dV^, we get the expression for the pre-symplectic structure in the 
{f3,4>,p^) parametrization as, 

^2 = ±-r^ — ■ .\n^ A d/3 + dp^ A d(P 
47rG7 sin(Ap)?7 

We can again restrict ourselves to > 0, which yields 

^2 = -r^ ^TV^^P't' d/3 + dp^ A d0 

AnG'yri sin(Ap) 

on the corresponding connected component. Regardless of the parametrization chosen, the 
pre-symplectic form Q has a degenerate direction corresponding to the Hamiltonian Vector 
field Xh^ of the Hamiltonian constraint H\. Physically the integral curves represent gauge 
direction along which points are physically indistinguishable. We can find the gradient of 
the constraint, 

VaH = 2r]^V sm{Xl3)^VaV + 2ri^\V^ sin(A/3) cos(A/3) V,/3 - 2p^V ^p^ 
and, using that = ^"''Vbif, and 

Q[a Qb] Q[a 



/ Q[a Qb\ Q[a Qb\ \ 



we get. 



X^lf = 167rG7r/Vsin(A/3)2 (^^^ - lQnGW>^V' siniX^) cosiXP) (^-^^ - Ap^ (J^^ 
as the restriction of the vector field to the constrained hypersurface with coordinates 



{V,P,(j)), using Eq. (A2). It is easy to check that this is indeed the degenerate direction 
of 0: Qi{Xh,-) = 0. 

We are now interested in understanding qualitatively the structure of the reduced phase 
space, for which we would need to take the quotient of T by the gauge directions generated 
by Xh- In this case, it turns out to be easier to perform a gauge fixing which, as we shall 
see, is well defined everywhere. Let us define C = — c, with c = constant. It is direct to 
see that it forms, with Hx a second class pair, 

{H, C}o = {H, 0}o = sm{X(3f - pj, 0}o = -{pj, 0}o = 2p^ 

provided 7^ 0, which is a condition we had asked before. The gauge fixed phase space 
r is therefore parametrized by {f3,p^) (or (V, (3)) with (p = constant, > or < 
and (3 G (— 1,0). The corresponding symplectic structures on each of this two different 
parametrizations are 

If we introduce the new variable (as was done in the main text), 

'X(3' 



X = In 



tan 



2 



then the symplectic form CI2 becomes 

Q2 = ±dp0 A do; 

which indicates that {x,p^) can be seen as conjugate variables in the reduced theory. 
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Appendix B: Some Integrals 

Here we list some integrals that are useful in the main part of the manuscript. 



x^'^+^e-^'dx = 0, with neZ 



X 



xo 



1 j 



d?7" y r] 



with n G Z 



Some error integrals. 



e ^ dx = ——erfc{z), with z & C 



f x^e ^'^dx = -e + ^erfc{z), with z e C 
where erfc{z) is the Complementary Error Function 

2 '•^ 



erfc{z) = 1 — erf{z) = 1 



e ^ dx 



which can be bounded using its asymptotic expansion [27] as 



\erfc{z)\ < 



-\z\^ 



-2 k/vr 



with |arg z\ < tt/A and z — )■ 00. When 2; = x is real then 



erfc{x) < 



Xi/tt 



1. Some Relations 



T]{k - (3f - 7]{k Tit- = -2r]R ik-aT 



irjt 
2r]~ 



u + a± 



trja 



trja 

u + a± =F ia 



+ 2ub± + b\ + 



a 
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Appendix C: Errors 



We show in explicit form the errors for the norm and volume and explain why the ap- 
proximations we take throughout the manuscript are completely justified. 



1. Norm Error 

The norm of the Gaussian states is given by 

/>oo 

Jo 

a 



- 2n+2 „oo 


V2ko 
u H 


J-^/2ko/cr 


a 



2n+l 



which for the n = case becomes, 



Ixlljhy = 2 



V2. 



V2ko 



/[u + y]e ^^du , with y = 



Then 



Ixlljhy = ^' 



■ue " dw + (j^ / ye " du 



+ CT yVTi - a y—erfc{y) . 



Using the form of y, we get 



Ixllphy 



27r/coO- + ye ^^-"'^^ - ak^J'^erfc 



'V2ko^ 



a 



The norm that we used in the manuscript was 

llxllphy ~ V27Tkoa. 
Then the error in this approximation is just 

,2 



Error(llxllphy) = ye ^ - ' ~ V 2 ^^"^^ I ~ 



(CI) 
(C2) 

(C3) 



(C4) 



(C5) 



(C6) 



This tells us that the error in the norm for n = is the error in the approximation of the 
function erfc{x) by the function ^r/^- If we use Eq. (B8) to bound the function erfc{y) 
then the error can be bounded by 



Error(||x 



""phy) < 2 ^ 



(C7) 



Since we have seen in the main text that semiclassicallity implies that we take /cq ^ cr > 
then this error is indeed very small. These errors were also studied numerically for n > 
and present similar behavior. We can then conclude that the approximation made is very 
good when computing the norm of the states. 
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2. Volume Error 



The coefficients in the expectation value for volume are, 



V+ 



K 



a 



(C8) 







/ 1 \ fcxi 



no 



1 n+1 



V2a V4 

with Mo = V^^o iia/ \/2. Then for n = are 



2u^ + AV2kou + Akl + e"" dw , (C9) 



± 



^ 1 A-oo 



2u^ + 4y2A;oM + 4A;^ + a^l g-^'/'^'dM . (CIO) 



Now we introduce other change of variables t = ^, and then the integral takes the form 



a 



2ah^ + 4:V2koat + Akl + e"* ^dt , (CI 1) 



with z = X ± = ^^^^^ ± The integral is 



/OO /* — 2; 

{2aV + AV2koat + Akl + e"*'dt - / (^aH"^ + 4^2/^0^^ + 4A;2 + j g-t^^^ 



4A;o + ^ - -^erfc{z) - 2V2(t- 



2A;o + ^ ) erfc{z) 



Then W takes the form 



K 



5a 



4/co + 



+ 



-2(7 _ 9 (T" 

-—^e 

'71 2y/nkQ 



erfc{z) - 2V2a 



a 



2A;o + ^ ) erfc{z) 



If we normalize using Eq. (|C5|) then 

K 



V- 



± 



4ko + 



+ 







re 



-^er/c(z) - 2V2a^-^ - (2^0 + ^ ) er/c(z) 



The values that we used in the manuscript was 



2a 



1 + 



+ 



(C12) 



Then the error in the integration is 



Error (Kj-) 



5a 



IT 



[za^ + 2V2a) + erfc{ 



+ 2A; + — 
2^kQ ° 2kQ 
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Using the triangular inequality 



I Error (14) I < — e^^'^e"'/'"' 
8a 



i\z\a^ + 2V2a) + \erfc{z)\ 



— ^— + 2ko + — 

2y/7ckQ 2/Co 



Using the relation (B7) (where |arg2;| = arctan ^ = arctan ^ < f because a <^ /ci 



|Error(l4)| < 



8^/7^a 



,±pQ^a2/2^2 |^|2 



a 



2ko 



2y/TikQ\z\ \z\ 2kQ\z\ 



< .J^2V2a^a^ko\z\e^''"^e^"'^^\-\-^\" 



1 1 

+ 



A^/2T{kl\z\'^a'^ 



+ 



+ 



2y/2a^kQ a'^akQ\z\ 
1 



2 I I 2^2 



|Error(\4)| < ^^^^^^kl ^±p»^aV2.^^-N^ 



As all quantities in the brackets are positive and less than one, we can bound the error as 

r\ 

2^2tx 

Given that V± can be too small or large then the quantity that really gives a measure of the 

ror(V 

|Error(V±)| 5a^cr^|z| 



error is the relative error H^EiEilk)! _ Ugj^g gq_ (|C12|), this becomes 



< 



/27r 



AkI 



< ^g-2fc2/a2g-a2/2a2 



27T 



The relation between the parentheses is less than one, then 



2 n1/2 



|Error(V±)| ^ 5a;^A;ocr^_2A:, 



^27r 



Since a ^ a, then e""^/^'^^ ^ 1. Furthermore, as fco ^ cr ^ a, then the negative exponential 
wins over the polynomial factor, which show that the error is small. These errors were also 
studied numerically for n > and present similar behavior. 

Appendix D: Tables 



In this appendix we summarize the main quantities that were found in the main text, for 
the generalized Gaussian states for n = 1 and compare them to the n = case. It should 
be noted, that for higher values of n, the structure of the terms is similar, with terms of the 



form 



ko 



and 



ko 
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I^llphy 



n=0 



2TTkoa 



n=l 



V± 



ae 



8ko 



ae 



±pa a2/2o-2 ( 3o-^+24fcgg-^+2a^cr^+16fc;^+8fcga^+a^ 



24fcocr^ 



Wo 



02 {kl + 



3o^ 
4 



02 



15(T'*+40fcg(T2+16A:4 



120-2 



w± 



a3e±2pae2aV-^ (A;2 + 3^2 ^ ^2) 



±2pa 2«Vo-^ / l5cr4+40fcgo-^+24a^o-^+16fcg+32a^A:g+16a 
"3^ e ( 12<T2+16fc2 



hkn 



1+ 



3<t2 



1+54 --^ 



15[T 

Idkf. 



3^2 
iff 



4 4fc2 



^12^2 



1+^ 



_9ct4 
16fcl 



16fe4 



4fc;^ 



ae 



.2/2.2 ^6^1 



(Ay)2 
(y>2 



,.2/^2 



,,2 2 



„a2/a2 



1 + 



3o^ 
4fc2 



1+^ 



2fen 



o 2 2 T 2 4 



1+- 



2^2 _ 




8fe; 
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